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ABSTRACT 

We report Giant Meterwave Radio Telescope (GMRT) continuum observations 
of six nearby normal galaxies at 333 MHz. The galaxies are observed with angular 
resolutions better than ~ 20" (corresponding to a linear scale of about 0.4 - 1 kpc). 
These observations are sensitive to all the angular scales of interest, since the resolution 
of the shortest baseline in GMRT is greater than the angular size of the galaxies. 
Further, for five of these galaxies we show that at 333 MHz, the mean thermal fraction 
is less than 5%. Using archival data at about 1 GHz, we estimate the mean thermal 
fraction to be about 10% at that frequency. We also find that the nonthermal spectral 
index is generally steeper in regions with low thermal fraction and/or located in the 
outer parts of the galaxy. In regions of high thermal fraction, the nonthermal spectral 
index is flatter, and has a narrow distribution peaking at —0.78 with a spread of 
0.16, putting stringent constraints on the physical conditions for generation, diffusion 
and energy losses of cosmic ray electrons at scales of ~ 1 kpc. 

Key words: galaxies: spiral — galaxies: individual (NGC 1097, NGC 3034, NGC 
4736, NGC 5055, NGC 5236, NGC 6946) — (ISM:) cosmic rays — (ISM:) Hii regions 
— techniques: interferometric — techniques: image processing 



1 INTRODUCTION 

The radio continuum emission from normal galaxies origi- 
nates from two emission processes: thermal free-free emis- 
sion from Hii regions predominantly seen at recent star 
formation sites, and the nonthermal synchrotron emission 
(hereafter referred to as nonthermal emission) due to accel- 
eration of cosmic ray electrons (CRe) in the ambient galactic 
magnetic field. A typical galactic magnetic field strength of 
few microgauss (1 - 15 /iG) and CRe in the energy range 1 
to 10 GeV give rise to radio nonthermal emission between 
0.1-10 GHz. However, only at frequencies below 1 GHz, does 
the thermal fraction of the emission reduce significantly, and 
hence low frequency studies are direct probes of nonthermal 
emission in galaxies. 

It is now abundantly clear that CRe in normal galaxies 
are accelerated in supernova remnant (SNR) shock fronts 
(with typical linear size < 1 pc) of Type II and Type lb su- 
pernovae produced by short lived (< 10^ yr) massive stars. 
Subsequently, the CRe diffuse away from their sites of origin 
and lose their energy primarily via synchr otron and inverse- 
Compton radiation in about 10* yrs (see I Condor] [l99'a for 
a review). In turbulent magnetic fields, such as those en- 
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countered in the spiral arms, the CRe diffusion speed is the 
Alfven speed (~ 100 km/s), and hence the CRe can expand 
to a volume of radius ~ 1 kpc without loosing much en- 
ergy, thus increasing the observed extent of the nonthermal 
emitting source. Note that the CRe would occupy a larger 
volume when observed at lower frequencies, since lower en- 
ergy electrons suffer lesser energy loss and hence can dif- 
fuse farther. The CRe generated at the shock front have 
steep electron energy spectra which results in a nonther- 
mal synchrotron spectrum, well represented by a power law, 
Si/.nt oc where Sv,-nt is the nonthermal radio flux den- 

sity at a frequency u and Ont is the nonthermal spectral 
index. The effects of CRe generation and propagation will 
cause Ont to vary from point to point in the galaxy. A t the 
acceleration sites one expects Ont ^ —0.5 to — 0.7 (|Belj 
ll978l : lBogdan fc VolkHTosl : iBiermann fc Strom|[T99l ). get- 
ting steeper due to various losses as the CRe diffuse away. 
To confirm the aforementioned scenarios a 

host o f studies have been don e in the past (see 

van der Kruit. Allen fc Rotd 
Beck fc Gravel 



e^ 



1977al 



jKlein et al. 


1982: 


Carilli et al. 


1993; 



Klein, Grave fc Wielebinski 
Niklas fc Beck 1997 : Berkhui jscn, Beck fc H ocrnea 
2003; Bec k et all boOsT [Murgia et all bood : 
iPaladino. Murgia fc Orrul |2009| ). However, due to the 
presence of thermal free-free emission, the Ont is con- 
taminated, making the measured value fiatter than the 
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Table 1. The sample galaxies observed with the GMRT at 333 MHz. 



Name 


RA 


Dec 


Morphological 


Angular 


i 


Distance 


Spatial 


Other data 








type 


size (D25)(') 


n 


(Mpc) 


resolution (pc/") 


available^ 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 


NGC 1097 


02h46ml9.0s 


-30dl6'30" 


SBbc 


9.3x6.3 


45 


14.5t 


70.3 


Ha, XlO/im, XlGOfim, 
A20cm 


NGC 3034* 


09h55m52.7s 


+69d40'46" 


lO 


11.2x4.3 


75 


3.97t 


21 


Ha, A20cm 


(M82) 


















NGC 4736 


12h50m53.0s 


+41d07'14" 


SAab 


11.2x9.1 


41 


4.661 


22.3 


Ha,A70/.tm, A160Atm, 


(M94) 
















A20cm 


NGC 5055 


13hl5m49.3s 


+42d01'45" 


SAbc 


12.6x7.2 


59 


9.2t 


44.6 


Ha, XlOfim, A160/im, 


(M63) 
















AlScm 


NGC 5236 


13h37m00.9s 


-29d51'56" 


SABc 


11.2x11 


24 


4.51^ 


22 


Ha, X70fim, A160^m 


{M83) 
















A20cm 


NGC 6946 


20h34m52.3s 


+60d09'14" 


SABcd 


11.5x9.8 


33 


6.8^ 


33 


Ha, X70fim, A160^m, 
A6cm, A20cm 



In column (5) D25 refers to the optical diameter measured at the 25 magnitude arcsec ^ contour from lde Vaucou leurs et al.| |l9 9lfl. 
Column (6) gives the inclinat i on an gle (i) defined such that 0° is face-on. Distances in column (7) are taken from: ^ iKarachentsev et al.l 
Karachentsev et "ahl ||2002| '1. ^ IKarachentsev. Sharina fc Huchtmeie j ll2000h and the NED Column (9) lists available ancillary 
data that have been used in this paper. ^ The wavelengths represents the waveband in which other data were available/used (see 
Section [2] for exact frequencies). * Thermal- nonthermal emission separation was not done for this galaxy (see text for details). 



true value. This is because the thermal emission has a 
spectral index of —0.1, while the nonthermal emission is 
significantly steeper (ctnt < — O.sj^- At frequencies below 
1 GHz, the nonthermal emission is expected to dominate. 
iNiklas. Klein fc Wielebinskil l| 19971 ) have shown for a sample 
of 74 galaxies, about 8% of the emission at 1 GHz is thermal 
in origin. At even lower frequencies, the contribution of 
thermal fraction reduces significantly, although there can be 
considerable local variation between, e.g, arms, interarms 
and giant Hll regions etc., so as to make the separation of 
the thermal emission important. Previous studies at low 
frequencies probed linear scales larger than 1 kpc, and hence 
any small scale structures were lost. These considerations 
emphasize the importance of high resolution, low frequency 
observations. 

In this paper we report 333 MHz interferometric ob- 
servation of six nearby galaxies, NGC 1097, NGC 3034, 
NGC 4736, NGC 5055, NGC 5236 and NGC 6946, using 
the Giant Meterwave Radio Telescope (GMRT). The galax- 
ies have angular size ~ 10' in the optical and the GMRT 
observations have a resolution better than 20", probing lin- 
ear scales of 0.4-1 kpc. The typical noise in the maps is 
~ 300 /^Jy/beam (except for NGC 3034), making them 10 
times more sensitive than previous stud ies at similar fre- 
quen c y of other nearby g a laxies ( Paladi no. Murgia fc Orr^ 
I2OO9I : [Heesen et aL| I2OO9I : ISukumar. Klein fc Gravd ll987D . 
All our maps have a dynamic range > 1000, which allows 
us to reliably determine the flux densities in regions near 
strong emission sites, mainly the nuclear region. We also 
use the technique developed bv lTabatabaei et al.l (|2007l ). us- 
ing archival Ha maps, for spatially separating thermal and 
nonthermal emission as discussed in Section[31 However, we 
emphasize that in using this robust technique, the resolution 
is limited by much coarser (40") far infrared maps (discussed 



1 Note that we use the parameter o for the total spectral index, 
where S^,tot « 1^", and Si/, tot is the total (thermal -I- nonthermal) 
flux density. 



in Section 4). In Section [2] we define the sample and discuss 
the observation and analysis procedure in Section 3. Section 
4 discusses the procedure to estimate Ofnt- Results on indi- 
vidual galaxies are presented in Section[5l We summarize the 
results in Section [B] The thermal nonthermal emission sepa- 
ration procedure using Ha and infrared data is discussed in 
Appendix [XI 



2 SAMPLE DEFINITION 

Our sample (see Table [l| includes six nearby spiral galaxies, 
with clearly demarcated arm and interarm regions as seen 
both in the optical and radio continuum. The optical size of 
these galaxies are typically 9'-12' and the GMRT synthe- 
sized beam (typically ~ 15") at 333 MHz can resolve the 
arm interarm regions. 

All the galaxies have 1.4 GHz integrated continuum 
flux density greater than 200 mjy (as measured in NVSS, 
ICondon et al.llT998l . data), which is necessary to ensure de- 
tection of low surface brightness diffuse emission across 
the galaxy. High resolution infrared and Ha data (see Ta- 
ble [1]) are also used. The infrared data were all taken 
from the Spitzer Infr ared Nearby Galaxy Survey (siNGS) 
iKennicutt et al.l |2003|) obtained usi ng Multiband Imaging 
Photometer (mips, Rieke et al.ll2004l ) on the Spitzer Space 
Telescope at A70 ^m and A160/xm, while the Ha images 
were obtained from the ancillary datasets of the SINGS and 
other public datasets, as described in Appendix [X] For five 
galaxies the data could be used for thermal nonthermal sep- 
aration, however in the case of NGC 3034, the A160/im mips 
image was not usable due to nonlinearity and streaking ef- 
fects. Finally, a higher frequency, near 1.4 GHz, archival 
continuum interferometric data/maps were used to find the 
spectral index with 333 MHz (see Table [1]). 

An interferometer with the shortest baseline Dmin, can 
detect all the flux from angular scales less than ~ 0.6A/Dmin, 
provided the UV-plane is densely sampled at the shortest 
spacings. The galaxies in our sample have optical angular 
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VLA UV coverage at 1465 MHz lor NGC 1097 



GMRT UV coverage at 330 MHz tor NGC 1097 
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Figure 1. UV coverage in the range of UV distance 0-220A for NGC 1097 at a declination of ~ —30.25°. The left plot is for the VLA 
in CD array and the right plot is for the GMRT. Sources of angular size above ~ 9' should be detected in this UV range. Note that the 
UV coverage for the VLA is sparse resulting in incomplete sampling of the UV plane, which leads to missing flux density problems above 
this angular scale. However, for the GMRT, several complete UV tracks are present. 



size less than ~ 12'. For face-on galaxies, the optical and 
radio continuum sizes are comparable, however, for highly 
inclined galaxies, the sizes may differ slightly. In GMRT 
^min ~100 meters, at 333 MHz allows us to detect all the 
flux from angular scales less than ~18'. Given that all the 
galaxies in our sample are less than ~12', the uv-plane is 
well sampled (see Figure [T]) and hence we do not expect any 
missing flux densitj^ This is an issue for the archival higher 
frequency data which is discussed below along with a short 
description of the sample galaxies. 

NGC 1097, in the optical is a spectacular barred spiral 
galaxy with prominent dust lanes in the bar. The central 
bar extends for about 20 kpc and then continues into the 
two optical spiral arms. The galaxy has an active Seyfert 
nucleus, and a circumnuclear ring of about 1.5 kpc diameter. 
Am ongst the major radio continu um studies of this galaxy, 
lOndrechen fc van der HulstI (|l983,) reported radio emission 
at 1465 MHz coincident with th e narrow dust lanes in the 
bar. Recently iBeck et al.l (|2005h obtained a spectral index 
map for this galaxy from observations at 4.8 and 8.4 MHz, 
which revealed relatively steep spectral index of —1 in the 
inner ridge of the bar, and becoming shallower to —0.7 as 
one moves to the outer ridge. 

To find the spectral index with GMRT 333 MHz, we 
analyzed the archival 1465 MHz data observed using CD 
array from the (project code: AW237). The map has 



^ The extent of missing flux density with our GMRT data was 
tested as follows: we took an existing archival 1.465 GHz map of 
NGC 6946, which was made by combining interferometric data 
from the Very Large Array (VLA) using C and D array jBeckl 
[2003). This galaxy is ~11' in size. The map was then sampled to 
the 333 MHz GMRT UV coverage and subsequently imaged. All 
the flux in the source (> 99%) could be recovered. 
^ The Very Large Array (VLA) is operated by the NRAO. The 
NRAO is a facility of the National Science Foundation operated 
under cooperative agreement by Associated Universities, Inc. 



an angular resolution of 40" x40", and the shortest measured 
baseline is ~150A, which can detect angular scales up to 
~14'. However, the uv coverage is sparse for angular scales 
above 8' (see Figure [l]). Since the size of the galaxy is about 
9', we expect missing flux density at large angular scales in 
these data. 

NGC 3034 (M82) is a prototypical, starburst, edge- 
on galaxy. The galajcy was initially classified as irregular, 
however recentlv iMavva, Carrasco fc Lunal (|2005l ). have re- 
ported the discovery of two symmetric spiral arms in the 
near infra-red. The galaxy is bright in the radio continuum, 
and has a synchrotron emitting halo extending up to a di- 
ameter of about 8 kpc. Based on a low frequency spectral 
index study of the whole g alaxy from 330 MHz to 4835 MHz, 
ISeaguist fc OdegardI (|l991j ) found a spectral index of —0.4 
in the nuclear region, steepening to —1 at a radius of 1 kpc 
due to transport of CRe from the disk by galactic winds. 

W e used the archival 1490 MHz VLA map (|Condonl 
1 19871 . obtained from the NED) at a resolution 60" x 60" to 
obtain the spectral index with 333 MHz. The 1490 MHz ob- 
servations were made with the VLA D configuration with 
shortest baseline length of ~ 170A, and are sensitive to an- 
gular sizes less than ~ 12'. This should be sufficient to detect 
the galaxy, which has an optical size (D25) of ~ 11'. How- 
ever, due to very bright emission from the core, the high 
resolution radio continuum map is dynamic range limited. 

NGC 4736 (M94) is the nearest large spiral galaxy, 
which has double ringed morphology. The inner ring is 
mainly made of recent star forming Hll regio ns, showing 
bremsstrahlung spectra (|Duric fc Dittmail[l9 88') and is dis- 
tinctly visible in the Ha. The outer ring however is a low 
surface brightness feature in the optical. The rings were sug- 
gested to be caused by inner and outer Lindblad resonances 
( Schommer fc Sullivanll 19761 ). The radio continuum study by 
Ide Bruvnl l | 19771 ) between 610 MHz and 1415 MHz, reported 
spectral index maps with the spectral index being ~ —0.5 



4 Basu et al. 



towards the center of the galaxy and steepening towards the 
outer parts of the galaxy. 

We used the archival WSRlQ 1374.5 MHz map_^taken 
from the Westerbork SINGS sample (,Braun et al. 2003) at a 
resolution of 19" x 12.5" after smoothing to 20" x 20" to 
obtain the spectral index. The observations was made using 
the "maxi-short" array configuration of the WSRT where 
the shortest east-west baseline was ~ 170A, 250A and 330A 
(36, 54 and 72 meters respectively), and are sensitive to 
angular scales less than ~ 12'. NGC 4736, having a angular 
size of ~ 6' in the radio, has good sampling of the uv-plane 
and thus this map is not affected by lack of short spacing 
measurements. 

NGC 5055 (M63) is a nearby flocculent spiral which 
has short multiple arms. Although this galaxy lacks orga- 
nized spiral arms in the optical, polarization observation at 
10.55 GHz shows regular, spiral magnet ic fields with radial 
component due to dynamo action (jKnapik et al-llioOOl ). Ra- 
dio continuu m spectral index maps be tween 610 MHz and 
1417 MHz bv lHummel fc Bosmal jl982l ) showed spectral in- 
dex of about —0.6 in the central regions of the galaxy, which 
steepens to —1 towards the outer parts. The galaxy appears 
featureless in the above study as well as in t he high frequency 
10.7 GHz map bv lKlein fc EmersonI (|l98ll ). 

The archival 1 8.5" x 12.5" 1696 MHz map from the 
Westerbork SINGS (|Braun et al.ll2007l ), after smoothing to 
20" X 20", was used for obtaining the spectral index with 333 
MHz. The "maxi-short" configuration of the WSRT at 1695 
MHz allows us to image angular scales up to 10'. However, 
the uv coverage is too sparse to detect flux from structures 
above ~ 7'. NGC 5055 has an angular extent >8', and thus 
we believe this map has missing flux density. 

NGC 5236 (M83) is a large barred spiral galaxy and 
is seen almost face-on. The central bar extends up to 3' 
(~ 4 kpc), showing clear dust lanes. The galaxy is bright 
in radio continuum and the central bar and the di s k is i n 
the 1465 and 4885 MHz observ ation of lOndrecheiil (|l985h . 
ISukumar. Klein &: Gravd (|l987l ) found the total and non- 
thermal integrated spectral index to be —0.75 and —0.8 re- 
spectively, obtained between 327, 1465 and 4750 MHz and 
a thermal fraction of 20% at 4750 MHz. 

The archival data at 1452 MHz VLA CD array (project 
code: AS325) were downloaded and analyzed to produce a 
26" X 14" map, which was used for obtaining spectral index 
with our 333 MHz observations. The 1452 MHz dataset had 
shortest baseline of ~ 150A which is sensitive to a maximum 
angular size of ~ 14'. The UV coverage is sparse for struc- 
tures above ~ 10'. NGC 5236 is about 11' in size and hence 
there may be a small amount of missing flux density in this 
observation. 

NGC 6946 is a large spiral galaxy with multiple arms. 
The galaxy has well separated arm and interarm regions 
and has been the subject of several studies in the radio 
continuum. There are several star forming regions in the 
galaxy and one of them shows the first likely detection of 
"anomalous" dust emission, due to spinning dust, outside of 



* The Westerbork Synthesis Radio Telescope (WSRT) is oper- 
ated by the Netherlands Foundation for Research in Astronomy 
(NFRA) with financial support from the Netherlands Organiza- 
tion for scientific research (NWO). 



the Milky Way (jMurphv et all bOld). Observations at 610 
MHz 1415 MHz and 5 GHz bv lvan der Kruit. Allen fc Rotj 
(Il977h revealed a bright radio disk of nonthermal origin. 
iKlein et all j 19821 ). based on observations at 10.7 GHz, es- 
timated a th ermal fracti on of 19 ± 10%. The spectral in- 
dex study bv lBeckl l|2007f ) between 1465 MHz and 8.5 GHz, 
clearly shows a flatter spectral index of —0.5 in the arms 
which steepens to about —1 in the interarm regions. 

For obtaining the spectral index with 333 MHz, we used 
an archival VLA 15" x 15" map at 1465 MHz. The VLA map 
was ma de by c ombining interferometric d ata from C c onfig- 
uration (|Beck| [2007') and D configuration l|Becklll99ll '). iBeckl 
(2007) noted the integrated flux density from the VLA map 
and single dish Effelsberg measurement had similar values 
and there was no indication of missing large angular scale 
structures. 

We emphasize that the GMRT 333 MHz images of the 
galaxies do not suffer from missing flux density, while some 
of the higher frequency maps (especially NGC 1097 and 
NGC 5055) may have missing flux density from large angular 
size structures. Thus, the true spectral index, particularly in 
the outer parts of the galaxies, is flatter than determined in 
this study. 



3 OBSERVATIONS AND ANALYSIS 

We conducted interferometric observations at 333 MHz for 
the six galaxies in our sample with the Giant Mctcrwave 
Radio Telescope (GMRT) near Pune, India (Swarup et al] 
Il99lf ). The data with 16 MHz bandwidth (corresponding to 
vb.f of 325 to 341 MHz, and a center ^'rf of 333 MHz), di- 
vided into 128 channels, were recorded following the usual 
protocol of observing flux and phase calibrator interlaced 
with observation on the source. At the beginning and end of 
the observing run, one of the flux calibrators 3C48, 3C286 or 
3C147 was observed for about ~ 10—20 minutes. At the be- 
ginning and end of the observing run, one of the flux calibra- 
tors 3C48, 3C286 or 3C147 was observed for about ~10-20 
minutes. Phase calibrators were observed every ~30— 35 
minutes for ~4— 5 minutes. Calibrators chosen for each of 
the observed sources are listed in Table (2] 

Data reduction was done using the Astronomical Im- 
age Processing System (aips) following standard procedure. 
After editing the data for strong radio frequency interfer- 
ence (RFI), standa rd flux and phase ca libration were applied 
to the source. The lBaars et al.l (|l977l ) absolute flux density 
scale was used to determine the flux densities of the flux 
calibrators and then applied them to the phase calibrator 
and the source. 

The task 'flgit' was used to remove low level RFI 
at 4 — 6(7 level across the frequency channels. Even lower 
level RFI was subsequently removed using the tasks 'tvflg', 
'wiper' and/or 'spflg'. The procedure was iteratively done 
by obtaining the gain solutions after every RFI removal 
stage. The flnal gain solution that was applied to the tar- 
get source when the closure errors were less than 1% on the 
phase calibrators. 

The calibrated data were then used to obtain the decon- 
volved images using the task 'imagr'. Prior to this, the data 
were condensed in frequency by vector averaging 6 adjacent 
channels of 125 KHz each (resulting in channel width of 750 
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Table 2. Observational summary of 333 MHz observations of our sample galaxies. Time spent on source excludes calibration overheads. 
Column (5) gives the synthesized beam achieved and column (6) gives the map rms (a"map). Column (7) gives the integrated flux of the 
source inside the icjmav contour. 



Name 


Obs. 


Time on 


Calibrators 


used 


Synthesized beam 




Integrated Flux 




date 


source (hrs) 






(" X ") 


{^Jy/beam) 


(Jy) 


(1) 


(2) 


(3) 


(4) 




(5) 


(6) 


(7) 


NGC 1097 


29-Oct-09 


6 


Flux Cal: 
Phase Cal: 


3C48 
3C286 
0116-208 
0409-179 


16x11 


300 


2.0±0.14 


NGC 3034 


31-Oct-09 


4.5 


Flux Cal: 


3C48 


22x15 


3000 


14±1 


(M82) 






Phase Cal: 


0834+555 
1459-1-716 








NGC 4736 


04-Jul-09 


7 


Flux Cal: 


3C286 


13x12 


2S0 


0.9±0.06 


(M94) 






Phase Cal: 


3C286 








NGC 5055 


02-NOV-09 


3 


Flux Cal: 


3C286 


17x10 


240 


2.3±0.2 


(M63) 






Phase Cal: 


30286 








NGC 5236 


31-Oct-09 


5 


Flux Cal: 


3C286 


16x12 


600 


7.4±0.5 


{M83) 






Phase Cal: 


1311-222 








NGC 6946 


29-Jul-09 


6.5 


Flux Cal: 
Phase Cal: 


30286 
30147 
1459-1-716 
2350-1-646 


12x11 


300 


4.3±0.24 



KHz) using the task 'splat'. This ensured that the band- 
width smearing was less than the size of one synthesized 
beam at 333 MHz. To account for wide field imaging with 
non-coplanar baselines, the technique of polyhe dron imag- 
ing was used in IMAGR (jCornwell &: Perlevll 19921 ). where the 
field of view is subdivided into a number of smaller fields 
(facets). For our purpose we used 7 x 7 = 49 facets covering 
the primary beam up to the half power beam width (HPBW 
~ 1.5°). Strong sources outside the HPBW, if present, were 
included in additional facets. 

Several rounds of phase only self-calibration were done 
iteratively, by choosing point sources such that, the flux 
density within one synthesized beam is more than 8a. Any 
stripe present in the map due to bad data, were removed us- 
ing a Fourier transform methocjfl. In the last iteration, one 
round of amplitude and phase self calibration was done. Final 
maps were made from full uv coverage and th e uv data wer e 
weighted using Briggs robust weighting of (|Briggslll995h . 
To CLEAN the extended diffus e emission from the galajcie s 
we used Sdi clean algorithm (|Steer. Dewdnev fc Itdll984 ). 
The final full resolution images obtained are shown in Fig- 
ure [2] and their corresponding synthesized beams are listed 
in column 6 of Table [21 

The integrated fiux density of the all the six galax- 
ies were obtained by integrating within the 4(Tmap (where 
o"map is the map noise) contour and is given in Tabled We 
compared the fiux density obtained for each of our sample 
galajcies with measurements done at other radio frequencies 
by other researchers, as given in Table [S] Figure [3] shows 
the integrated broadband spectrum for our sample galaxies. 
Our flux density measurements are in good agreement with 
the interpolated fiux densities from higher and lower radio 
frequency observations reported in earlier studies. 

The uncertainties in the estimated flux density (S'source) 
depend on the rms noise in the map as well as on errors as- 
sociated with uncalibrate d system tempera ture (Tsys) varia- 
tions, which is about 5% (|Rov fc Radl2004l 'l at 333 MHz for 



^ see |http://www.ncra.tifT.res.in/~aritra/stripe-removal.pdf| 



the GMRT. The fiux density error (5Ssource) can be calcu- 
lated as. 



^•S'source — \ / ( T' ^ '5'source | ^" f^map 

The fiux density also has an systematic error of about 8% 
associated with the absolute fiux scale error at 333 MHz. 



4 NONTHERMAL SPECTRAL INDEX (ont) 

The nonthermal spectral index ant, is the main parameter 
of interest because it is used to model CRe generation and 
propagation. However, it is not possible to measure it di- 
rectly. The quantity which is measured is the total spectral 
index a which is contaminated by thermal free-free emis- 
sion. This contamination is significant for spatially resolved 
regions of the galaxy where the thermal fractiorjj (hence- 
forth, /th) is high. In this section, we assess the effect of 
fth on Ont obtained between 333 MHz and near 1.4 GHz 
in spatially resolved regions of the galaxy. To do this, we 
use the method of thermal-non thermal separation devel- 
oped by iTabatabaei et al.l |20o3) and apply it to five of our 
galaxies (except NGC 3034). The method uses dust extinc- 
tion corrected Ha maps to get a template for the free-free 
emission across the galaxy, and extrapolates it to the de- 
sired radio frequency to obtain spatially resolved thermal 
radio emission maps. The extinction correction is achieved 
by estimating the dust temperature and optical depth us- 
ing the far infra-red maps at A70/im and A160/im (see Ap- 
pendix |X] for details). The far infrared XlGO^m map has 
an angular resolution of 40" x 40", and hence the final 
nonthermal spectral index maps have this angular resolu- 
tion, which is significantly coarser than the total spectral 
index (hereafter a) maps. The usefulness of this method lies 

^ Thermal fraction is defined as: /th = S'^,th/'S'i/,tot, where, S„ tj, 
is the flux density of the thermal emission. In the text we express 
/th as percentages. 
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Figure 2. The contour maps arc GMRT 333 MHz observations of the sample galaxies which are overlaid on the optical DSS image in 
greyscale. The contour levels starts from 4cr, increasing in multiples of \/2. The grey contours shows (—2, —3, —4) X a. 
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Figure 3. Integrated flux density as a function of frequency. The black circles represent the integrated flux densities without the thermal- 
nonthermal emission separation (see Table [Sjl, while the red triangles represent the nonthermal flux densities of the galaxies. The blue 
square for NGC 5055, is the flux density measured within 6.7' X 3.2' (see Section|4]l. 



in the fact that it can b e used in spatiall y resol ved parts 
of the galaxy. Further, iTabatabaei et aP |20o3) demon- 
strate the robustness of their method over several other 
existing thermal-nont hermal separation techniques (see ref- 
erences therein and Broadbent. Osborne fc HaslamI [l989l : 
iGioia. Gregorini fc Kleinlll982l : iHummel fc BosmalllOSd V 

The mean fth for all our sample galaxies was found to 
be less than 5% at 333 MHz, although the value can go up 
to 10% in some specific bright Hll regions and in the spiral 
arms. Applying the same method to the available higher 
frequency data at or close to 1.4 GHz, we found the mean 
fth to be less than 12%. The distribution of the fth at 333 
MHz and near 1.4 GHz for five galaxies is shown in Figure lAl] 
of the appendix. 

The nonthermal 333 MHz and higher frequency maps 



were further used to obtain Qnt variation for the five galaxies. 
The left panels of Figure 31 show the a maps and the right 
panels show the ant maps. Figure [5] shows the variation of 
a and ant with radius for the galaxies. Here, the spectral 
index was estimated by azimuthally averaging over annuli 
of one beam width. The radial profiles clearly reveal that 
the nonthermal spectral index is nominally steeper than the 
uncorrected values, and further steepens towards the outer 
parts of the galaxy. However, we note that towards the outer 
parts of the galaxies, the steepening of both the total and 
nonthermal spectral indices may be caused due to "missing 
flux density" problems. 

We are now in a position to quantify the importance 
of thermal-nonthermal separation, for robust estimates of 
ant between 333 MHz and near 1 GHz. The uncertainty in 
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Figure 4. Left column: The total spectral index (o) maps between GMRT 333 MHz and near 1 GHz. The top row is NGC 1097, middle 
row is NGC 3034 and the bottom row is NGC 4736. The a maps for NGC 1097, NGC 3034 and NGC 4736 have a resolution of 40" X 40", 
60" X 60" and 20" X 20" respectively. Right column: The nonthermal spectral index (ont) between 333 MHz and near 1 GHz at resolution 
of 40". Overlaid are the 333 MHz contours. Contour levels are indicated below each figure. 
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Figure 4. contd... Left column: The total spectral index (a) maps between GMRT 333 MHz and near 1 GHz. The top row is NGC 5055, 
middle row is NGC 5236 and the bottom row is NGC 6946. The a maps for NGC 5055, NGC 5236 and NGC 6946 have a resolution of 
20" X 20", 26" X 14" and 15" X 15" respectively. Right column: The nonthermal spectral index (ont) between 333 MHz and near 1 GHz 
at resolution of 40". Overlaid are the 333 MHz contours. Contour levels are indicated below each figure. 
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Figure 5. Azimuthally averaged spectral index determined within annuli of one beam width. The black circles are the total spectral 
index (o), while the red triangles are the nonthermal spectral index (ctnt)- The top x-axis shows the angular distance from the center in 
arcmin. For NGC 3034, the averaging was done for half the synthesized beam and therefore adjacent points are not independent. 
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Figure 6. The bottom panel of the plot shows the ratio a/ant as 
a function of /t^ at 1465 MHz for NGC 6946 obtained at spatially 
resolved regions of the galaxy. The black circles corresponds to 
o/ont between 333 MHz and 1465 MHz, while the red triangles 
shows between 1465 MHz and 4860 MHz. The top panel shows 
the flux density normalized to the peak value at 333 MHz (in 
black) and 1465 MHz (in red) respectively. Note that generally 
higher flux density regions have high thermal fractions. 



Ont can be estimated by propagating the various sources 
of errors in obtaining the nonthermal maps. Three primary 
sources of error affect the ant measurement, namely: the rms 
in the two radio maps (o-^.map), the error associated with the 



uncalibrated system temperature of the instrument (T^,Tsys, 
and the error in thermal fraction, ctt. f^^Tsys is about 5% at 
333 MHz for the GMRT and assumed about 2% for VLA and 
WSRT near or at 1.4 GHz maps. The err or <jt is estimated 
to be about 10% at these frequencies by iTabatabaei et abl 
(|2007h (also see Appendix A). Incorporating these errors, 
we have computed the ratio a/ont in regions by dividing 
the galaxy into rectangular grids of size of approximately 
40". We find that for thermal fraction less that 5% at 333 
MHz and less than about 10% near 1.4 GHz, a-nt ^ a. This 
effect is illustrated for NGC 6946 in the bottom panel of 
Figure 6, where the a /a-nt is plotted (in black) as a function 
of thermal fraction near 1.4 GHz. Within the error bars the 
ratio is consistent with unity for thermal fraction less than 
10%. For higher thermal fraction (predominantly in the arms 
or Hll regions) the ratio drops below unity, and correction 
to a, to obtain a-nt, could be as high as 20%. The thermal- 
nonthermal separation is even more important for ant cal- 
culated between higher frequency pairs like 1.4 GHz and 
4.8 GHz. The ratio a/ ant for NGC 6946 obtained between 
1465 MHz and 4860 M Hz (using archival map downloaded 
from NED, lBeckll2007l ) is plotted in red in Figure 6, where 
we clearly see that the ratio is significantly below unity for 
the whole range of thermal fractions. Note that since there 
is also a systematic uncertainty associated with the abso- 
lute calibration scale for calibrating radio flux densities, the 
absolute value of the spectral index can vary by about 10%. 



5 DISCUSSION ON INDIVIDUAL GALAXIES 

NGC 1097: Our 333 MHz map (see Figure [l]) has an an- 
gular resolution of 16"xll", which corresponds to a spatial 
resolution of ^ 1.1 kpc, if the galaxy is at a distance of 14 
Mpc. The only other low frequency radio continuum study 
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Table 3. Multifrequency integrated flux density for our sample 
galaxies. 



Source 


Frequency 


Flux density 


Ref. 




(GHz) 


(Jy) 




NGC 1097 


0.333 


2.1±0.14 


This paper 




1.365 


0.350±0.025 


1 




1.465 


0.42±0.02 


VLA CD array 




1.665 


0.262±0.030 


1 




4.850 


0.142±0.012 


1 




8.450 


0.094±0.018 


1 


NGC 3034 


0.038 


23±3 


2 












0.333 


14±1 


This paper 




0.750 


10.7±0.5 


2 




1.415 


8.0±0.4 


3 




2.695 


5.7±0.3 


2 




5.0 


3.9±0.2 


2 




10.7 


2.25±0.06 


4 




14.7 


1.79±0.04 


4 


NGC 4736 


0.333 


0.9±0.06 


This paper 












0.610 


0.435±0.026 


7 




1.365 


0.320±0.01 


6 




1.415 


0.276±0.02 


7 




4.995 


O.llliO.Ol 


7 




10.7 


0.096±0.018 


8 




15 


0.017±- 


9 


NGC 5055 


0.333 


2.3±0.13 


This paper 




0.408 


1.16±0.022 


5 




0.610 


0.83±0.04 


10 




0.750 


O.S8±- 


11 




1.417 


0.409±0.02 


10 




1.696 


0.285±0.02 


WSRT 




2.295 


0.26±0.02 


12 




3.7 


0.254±- 


13 




10.7 


0.117±0.026 


14 


NGC5236 


0.327 


6.86±0.62 


15 




0.333 


7.3±0.4 


This paper 




1.452 


2.3±0.15 


VLA CD array 




1.465 


2.2±0.11 


15 




2.7 


1.24±- 


16 




4.75 


0.97±0.13 


15 


NGC 6946 


0.158 


7.54;-- 


17 




0.178 


5.88±0.44 


18 




0.333 


4.e±0.24 


This paper 




0.408 


3.95 


5 




0.610 


2.34 


19 




1.415 


1.54±0.10 


19 




1.465 


1.41±0.08 


VLA C + D array 




2.695 


0.83±0.03 


20 




4.86 


0.54±0.05 


VLA+Effelsberg 




8.50 


0.45±0.05 


21 




10.70 


0.359±0.042 


8 



IBocIc et al.' J2002f), ^Kellcrmann et al.' flOSgl), ^Hummc? I'lgSO"), 
^ Klein, Wiolebinski & Morsi (1988), ° Gioia & Gresorini 11980), 
fe raun et al. 12007), "'do Bruyn (1977), ^'Klein & Emerson (1981), 
^Nag ar. Falcke & Wilson (2005), l'-" Hummel & Bosma (1982), ^Me .JonR 
<196a l, l^Arp (1973), Kuril'chik, Onishchenko & Turyvskii ( 1967i ), 
^%ukumar, Klein & Grave (1987), ^'"Wright et al. (1990j), 
^^.Brown fc Hazar d 11961). Caswell & Wills (1967), 
^ Hvan der Kruit. Allen & Rots ( 1977), ^ Ide Jongl jl967l ), 
^ iKuril'ch'ik et al. ( 1970l). 



of th is ga lcixy was done by Ondrechen fc van der HulstI 
l|l983h and lOndrechen. van der Hulst fc Hummell (| 19891 ) at 
1465 MHz. They reported strong radio emission coinci- 
dent with the narrow dust lanes in the bar. The central 
bar region is also clearly visible in our 333 MHz map, 
and is coincident with dust emission inferred from the 
far infra-red A24^m images from Spitzer, as seen in Fig- 
ure 7. The bar pattern rotates in the clockwi se direction 
(see I Ondrechen. van der Hulst fc HummelllT989l ). and as ex- 
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Figure 7. The bar in the galaxy NGC 1097. In shades of orange 
is the 24/im MIPS image (resolution ~6"), overlaid is the high 
resolution 333 MHz map made using baselines more than 1 kiloA 
with a resolution 14.5" X 10". The contours start from 1 mjy, and 
subsequent contours are in multiples of V2. The highest intensity 
contour is at 150 mjy/beam. 



pected prominent shock fronts are visible on the leading side 
of the bar with respect to the sense of rotation. The fth in 
the bar region at 333 MHz is estimated to be <1%, hence 
the radio emission at this frequency is entirely nonthermal in 
origin. These observations fit well with the explanation that 
the radio emission from the bar results due to shock com- 
pression of gas, dust and magnetic field, as was concluded by 
lOndrechen fc van der HulstI (jl983l ). The total emission from 
the bar and the central component (the region shown in Fig- 
ure 7) is 835±57 mJy, out of which 750±52 mJy arises from 
the central co mpact region (25"x2l"). For a sim ilar sized 
central region lOndrechen fc vaxi der HulstI l)l983l l quote a 
flux density of 260 mJy at 1465 MHz, based on which we get 
a a of ~ —0.75. The a^t within a 40" region based on our 
analysis is —0.85 ± 0.1. This Qnt close to the nuclear region 
is steeper when compared with other galaxies in our sam- 
ple. This indicates possibility of radiative cooling by inverse 
Compton/synchrotron processes. 

The southeastern (SB) side of the bar has a flux den- 
sity of 50 ± 4 mJy and the northwestern (NW) side has a 
flux density of 30 ±2 mJy, measured within the 2 mJy closed 
contours. The red lines in Figure[7]show the SE bar and NW 
bar. In the optical, the bar region continues into two promi- 
nent arms however the radio emission is highly diffuse. We 
flnd the thermal fraction of the interarm region at 333 MHz 
and 1465 MHz to be 0.2% and 2.5% respectively whereas in 
the spiral arm regions it is 1% and 4.5% respectively. In the 
arm regions, the (ont) is —1.4 ± 0.09 while it gets signifi- 
cantly steeper (—1.8 ± 0.12) in the interarm regions which 
may be caused due to missing flux density problems. The 
radial proflle of the a and ant in Figure (5] gets steep to- 
wards the outer parts of the galaxy possibly due to missing 
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flux density. We emphasize that the overall thermal fraction 
in this galaxy is very low, with (/th) of ^^0.6% at 333 MHz 
and <4.8% at 1465 MHz (see Figure [AHi. 

At 333 MHz , we detect the radio galaxy discovered by 
iBeck et al.1 (j2005l l. located ~90" from the nucleus of NGC 
1097 towards the southwest at a PA of -110° (see Figure 7). 
The radio galaxy has a total flux density of 9.5 ± 0.5 mjy, 
with 4.8 ± 0.25 mJy in the northern component (RA = 
02h46ml2.7s, dec = -30dl7'01" (J2000)) and 3.8 ± 0.2 
mJy in the southern component (RA = 02h46ml2.5s, dec = 
-30dl7'26.6" (J2000)). 

NGC 3034 (M82): The 333 MHz map (see Figure [g]) 
of this galaxy has an angular resolution of 22" x 15", cor- 
responding to a spatial resolution of about 0.5 kpc for a 
distance of 3.9 Mpc. We flnd the integrated flux density at 
333 MHz to be 14 Jy which is consistent with flux densi- 
ties measured at other frequencies. The rms of the image, 
which is 3 mjy/beam, is dynamic range limited since the 
peak flux density in the galaxy is 3 Jy/beam. Morphologi- 
cally the galaxy is featureless at our resolution, although the 
northern extension (RA = 09h55m56s, dec = -|-69d4l'48") 
is coincident with the base of the Ha/optical horn (see Fig- 
ure [2|, which is likely due to synchrotron emitting particles 
convected outwards by the nuclear wind, tracing the Ha 
horns. The galajcy was observed be tween 330 MHz to 4835 
MHz bv lSeaguist fc OdegardI (|l99ll ). and a similar horn was 
seen at ~1515 MHz. Their study also revealed a nonthermal 
radio halo surrounding the galaxy, however our observations 
do not show this due to dynamic range limitations. 

For this galaxy, the thermal-nonthermal separation was 
not possible, but we have computed the a map and radial 
profile as shown in Figure H and [S] The 333 MHz GMRT 
map was convolved to the resolution of VLA map at 60". 
This is a good approximation for Ont index variation across 
the galaxy since the thermal fraction of this galaxy has been 
estimated to be 15% at 32 GHz fKlein. Wielebinski fc MorsH 
1 19881 ). corresponding to ~ 3% and 2% at 1.4 GHz and 333 
MHz respectively. The spectral index is seen to be very flat 
with a = —0.35 ± 0.03 towards the center and steepens to 
— 1 ± 0.2 at a distance of ~2 kpc from the center. Similar 
radial variation of the spe ctral index was also reported by 
ISeaauist fc OdegardI (|l99ll ). 

NGC 4736 (M94): The 333 MHz map (see Figure [H) 
of this galaxy has an angular resolution of 13" x 12", corre- 
sponding to an spatial scale of about 0.3 kpc for a distance 
of 4.7 Mpc. The other low freque ncy study of this galaxy ex- 
ist at 610 MHz (|de BruvnII 19771 ) . observed with much lower 
angular resolution of 56" x 85". At our angular resolution we 
clearly detect emission from the inner star forming circum- 
nuclear ring of the galaxy as shown in Figure 8 overlaid on 
the Ha map. The ring is located at about 50" (major axis), 
corresponding to 1.1 kpc, from the center, and is ~19" (420 
pc) wide (up to the half power) having a flux density of 
240±13 mJy. The estimated thermal fraction of the ring at 
333 MHz and 1374.5 MHz is 5 ± 1% and 10 ± 4% respec- 
tively, and hence the ring at 333 MHz is largely nonthermal 
in origin. The ring is seen to be coincident with the Ha ring, 
which is shown in Figure[S]after convolving to the same radio 
resolution. 

Radio emission from the ring is also see n in the higher 
i reque ncy 8.46 GHz observations made by IChvzv fc Butal 
120081 ) ■ however at these frequencies the emission is a mix- 
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Figure 8. The star forming inner circumnuclear ring in NGC 
4736. The contours are the high resolution 333 MHz GMRT map 
with synthesized beam of 10.8" x 9.9" and the grey scale shows the 
Ha image (in erg s~^ cm~^) from the Jac obi Kapteyn Telescope 
with the 6570A filter jKnapen et"alll2004) . smoothed to the same 
resolution as the 333 MHz map. The lowest contour is at 3 mJy 
increasing in steps of 1 mJy. 



ture of thermal and nonthermal components. The ring has 
relatively high thermal fraction which leads to a noticeable 
change in spectral index (see Figure[5]) : (a) is —0.45 ±0.02, 
while the (ant) is —0.6 ± 0.1. The ring has massive on going 
star formation dSmith et al.| [l991) and ant value indicates 
that both CRe generation and escape happens on timescales 
shorter than the radiative timescale. The compact nucleus 
lies within the central 0.5 kpc (25"x22") with integrated 
flux density of 110±6 mJy and a peak flux of 45 m Jy/beam 
and ant = —0.6 ± 0.05. Outside the ring extending to the 
edge of the galaxy there is no thermal emission observed 
leading to ant ^ a, 

NGC 5055 (M63): The 333 MHz observation (see 
Figure (21) has an angular resolution of 17"xl0" correspond- 
ing to a spatial scale of 0.7 kpc at a distance of 9 Mpc. 
As is evident from the flgure, the extent of the galaxy is 
larger at 333 MHz than the optical. In an earlier study by 
iHummel fc Bosmal J 19821 ) at 610 MHz and 1410 MHz at 
the much lower angular resolution of 58"x87", the galaxy 
appears featureless which is also the case in our higher reso- 
lution 333 M Hz map. However, 1 0.55 GHz polarization ob- 
servations bv lKnapik eTell (|2000h show regular spiral struc- 
ture. The overall radio emission at 333 MHz at a level of Aa 
is ~8.6' x5.3' in extent, which has an integrated flux density 
of 2.3±0.13 Jy. This value is higher than what one would 
obtain using a spect ral index of —0.78 estim ated between 
610 and 1417 MHz ("Hummel fc Bosmal [l98^ ) . The reason 
for this is lHummel fc Bosma (1982) estimate the integrated 
flux density within a region 6.7' x3. 2' at 610 MHz, which 
is smaller than our detected size at 333 MHz. However, if 
a similar size region is considered, we obtain a flux density 
of 1.4±0.1 Jy which is consistent with a spectral index of 
—0.78 as seen in Figure [3|d). NGC 5055 is highly inclined, 
with inclination angle, i — 59°, and the extended low fre- 
quency radio emission can result due to escape of low energy 
(< 2 GeV) CRe from the disc which travel larger distances 
4.5 kpc) into the galactic halo. 
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Central region of NGC 5236 




Right Ascension (J2000) 

Figure 9. The central region in NGC 5236. Greyscale is the DSS 
image at a resolution of 1", overlaid are the 333 MHz contours at 
resolution 13" X 10". The marked regions are as per Ondrechen 
(1985). The flux densities of the sources are as follows; A : 12.3 it 
0.7 mjy, D : 11.5 ± 0.7 mjy, F : 13.5 ± 0.8 mjy, G : 55 ± 3 mjy, 
H : 10.2 ± 0.6 mJy, I : 9.2 ± 0.5, X : 20.0 ± 1.2 mjy 

The mean thermal fraction (/th) at 1696 MHz (A18cm) 
and 333 MHz is 12% and 1.3% respectively. The Ont in 
the central region is ~ —1 which steepens to —1.7 towards 
the outer parts of the galaxjQ ( see Figure [5l) . This result 
supports the conclusion drawn bv lHummel fc Bosmal (|l982l ) 
that the spectral index steepening is mainly due to energy 
losses of CRe and decrease in number of relativistic electrons 
with increasing galactic radius. 

NGC 5236 (M83): The 333 MHz observation (see 
Figure [2]) has an angular resolution of 16"xl2", which 
correspond to a linear scale of about 0.4 kpc at a dis- 
tance of 4.2 Mpc. An early study of this galaxy by 
ISukumar. Klein fc Grav3 (| 19871) at 327 MHz using the Ooty 
Synthesis Radio Telescope (OSRT) quotes integrated flux 
density of 6.86±0.62 Jy at a resolution of 53" x 33", while 
our GMRT observations has a slightly higher integrated flux 
den sity of 7.4±0 4 Jy 

lOndrechenI (| 19851 ) has reported several point sources in 
this galaxy. To verify the presence of these sources at 333 
MHz, we made a high resolution i mage of 13" x 10" a s shown 
in Figure [51 Comparing this with lOndrechenI (|l985l ) we de- 

Note that iHummel fc Bosmal lll982l ') using observations at 610 
and 1417 MHz quote a spectral index variation of —0.6 to — 1 
between the central and outer parts of the galaxy, while we find 
much steeper values. The WSRT 1696 MHz map has a flux den- 
sity of 284 mjy within the same region of 6.7'x3.2', which when 
extrapolated using a spectral index of —0.78, has a much lower 
value (see Figure |3jd) ) . This may be due to missing flux den- 
sity issues as discussed in Section [2] Thus the true spectral index 
should be flatter than what is observed by us. 



tect the central bar (shown in red lines) and the Hll regions 
named as A, F, H and I, and nonthermal polarized point 
sources D, G and X. Note that source G is identified as a 
background galaxy by ,Maddox et al.l (|2006l ) . The region E 
was identified as a shocked region in the arm, and has a fiux 
density of 33 ± 2 mJy, while the point source Q is a Hll re- 
gion and has fiux density of 5 mjy/beam. However, we were 
unable to detect point sources B and C. The fiux densities of 
the point sources are given in Figure In addition to this 
we see an extended radio arm, which runs along the faint 
narrow dust lane visible in the infrared A70/^m map. 

The galaxy has a low mean thermal fraction of about 
3% at 333 MHz and about 7% at 1452 MHz, and hence 
Qnt ~ a as seen in the radial spectral index profile of the 
galaxy in Figure [5] In the arms, at 333 MHz, the thermal 
fraction is estimated to be 8 — 10%, where a ranges from 
—0.3 ± 0.04 to —0.5 ± 0.05 which gets modified to ct-nt as 
—0.4 ± 0.05 to —0.65 ± 0.05. Qnt along the arms is seen to 
change significantly from —0.4 within the central 3 kpc to 
—0.65 beyond 3 kpc. Thermal fraction at 333 MHz in the 
interarm regions is 1 — 4% and beyond the central 4.5 kpc is 
less than 1%, indicating purely nonthermal emission. In the 
interarm regions, a-nt values lie in the range —0.7 ± 0.06 to 
-1.2 ±0.09. 

NGC 6946: The 333 MHz observation (see Figure 
has an angular resolution of 12"xll" which corresponds to 
a spatial resolution of 0.4 kpc at a galaxy distance of 6.8 
Mpc. The integrated fiux density at 333 MHz is 4.5±0.24 
Jy. Using our 333 MHz and archival 1465 MHz VLA maps, 
we estimate the mean thermal fraction to be about 4.3% 
and 11% respectively. 

Spectral index study bv lKlein et al. I l|l982l ) between 610 
MHz and 10.7 GHz showed variation from about —0.55 near 
the center to about —0.9 towards the edge. The spectral in- 
dex between GMRT 333 MHz map and VLA 1465 MHz 
map shows a similar trend, with spectral inde x of — .45 at 
the center and —1 towards the edge. Recently, iBeckl (|2007| ) 
studied the spectral index between 4860 MHz and 1465 MHz 
which showed flatter spectral index in the arm than the in- 
terarm due to higher thermal fraction. In our studies, ther- 
mal fraction is higher in the arm (5-10%) than in the in- 
terarm (1-5%). In the arm, a lies in the range —0.55 to 
—0.65, while the Qnt is steeper —0.7 to —0.85. In the in- 
terarm regions the average spectral index is steeper with 
(a) = —0.9 ± 0.06, while the Qnt does not change signifl- 
cantly, {ont} = -0.95 ±0.11. 

The 333 MHz map prominently shows the g iant H ll re- 
gions identified bv lvan der Kruit. Allen fc Rot^ |l973). The 
Hii region at RA = 20h35m06s, dec = ±60dl0'5 6" (J2000) 
(region C in Ivan der Kruit. Allen fc Rots! l| 19771 )). is about 
35" (1.2 kpc) in size, and can be seen in the 40" low res- 
olution nonthermal maps. This region is a suitable candi- 
date to investigate the CRe injection near high star forming 
sites. We find the thermal fraction to be 10 ± 3% at 333 
MHz and 15 ± 4% at 1465 MHz. The Qnt in this region is 
— 0.6±0.08, whereas before the separation the spectral index 
is -0.4 ±0.05. 

The central 35" region (at RA = 20h34m52s, dec = 
±60d09'l4" (J2000)) has a flux density of 226±13 mJy at 
333 MHz while at 1465 MHz the same region has a fiux den- 
sity of 117±6 mJy. This corresponds to an a of —0.44 ±0.05 
at a resolution of 15" x 15". However, after the separation. 
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similar region has a Ont of —0.54 ± 0.1 (see Figure |4]). In 
the nuclear region about 15% of the total emission is due to 
thermal emission at 333 MHz, higher than the rest of the 
galaxy, where the thermal emission lies within 10%. 



6 SUMMARY 

The 333 MHz GMRT observations of six nearby galaxies 
(greater than 10' in size) presented in this paper are by far 
the highest angular resolution maps for these galax;ies at 
this frequency which are sensitive to both small (~ 0.5 kpc) 
and large (>10 kpc) scale structures. The hybrid array of 
GMRT enables us to recover emission from both the diffuse 
and small scale emission in these sources, which we have ad- 
equately verified by comparing our estimated integrated flux 
densities with those available in the literature. Our observa- 
tions have comparable resolution as that attained at higher 
frequencies, and probe linear scales of about 0.4-1 kpc. 

At 333 MHz the galaxies appear smoother than at 
higher frequencies, where spiral arms are often discernable. 
The thermal emission (which is primarily associated with 
ionized structures observed in Ha maps) for all the galaxies 
was found to be less than 5% of the total emission, and hence 
at these frequencies, the emission is essentially nonthermal. 
The smooth appearance is consistent with the conjecture 
that nonthermal emission at 333 MHz primarily results from 
a population of old (> 10^ yrs) CRe, which diffuses away 
from their formation sites filling up a volume of radius > 1 
kpc, without loosing much energy. This scale is larger than 
the width of the gaseous spiral arms and comparable to the 
width of the inter-arm region (see e.g. Condon 1992), hence 
reducing any contrast in intensity across the galaxy. The 
prominent radio structures observed at higher frequencies 
(> 1 GHz) are visible solely due to the increased thermal 
fraction. 

We have determined robust estimates of spatially re- 
solved Qnt (see Section ^ for five galaxies in our sample, 
with a resolution of 40", corresponding to a linear scale of 
1-3 kpc, which also corresponds to the diffusion scale of 
the CRe. The observed Qnt ranges from —0.3 to —1.8 in 
various regions of the galaxies. Generally CRe are thought 
to be a ccelerated at SNR where both the theoretical pre- 
diction (|Belllll97^ ) and ob servation based on Galactic SNR 
suggest mean ant ~ -0.5 (|Greenlll99^ : lKothes et al.ll2006f ). 
The SNR are short lived (~ 10^ yrs), while the CRe dif- 
fuse away from their acceleration sites losing their energy 
through several physical mechanisms for typically 10* yrs, 
and hence their initial energy spectrum gets distorted lead- 
ing to a change in Qnt . If the en ergy losses are due to ioniza- 
tion, Qnt flattens (jLongair 201l'). Synchrotron radiation and 
inverse Compton scattering loss leads to a steepening of the 
spectra, while adiabatic cooling and bremsstrahlung keep 
the spectra unchanged. To assess the behaviour of Qnt in 
spatially resolved parts of the galaxy we plot Qnt distribution 
(as shown in Figure fTO)) for two regimes of thermal fraction: 



fit > 10% (filled grey histogram) and fit < 10% (un- 
filled histogram), using data from 4 galaxie^. To do this. 



° We did not use NGC 5055, since we think that there is a missing 
flux density problem at 1690 MHz and hence the derived a and 
Ont might be steeper than the actual value. 



{O^^^ > 10%) 



-1 — I — I — I- 




(f,/™' < 10%) 



-0.5 



-1.5 



Nonthermal Spectral Index (a^^j) 



Figure 10. The above plot shows two histograms of nonther- 
mal spectral index for 4 galaxies for two range of thermal frac- 
tions. The filled grey histogram corresponds to thermal fraction 
/th*°"^ > 10% and the unfilled histogram is for flt'^^^ < 10%. 



we have computed Qnt by dividing the galaxy in rectan- 
gular grids of 40". The histograms are binned with step 
size 0.05 in Qnt. The high thermal fraction corresponds to 
bright regions in Hq which traces gas ionized by massive 
OB stars in star forming regions. These regions host the 
CRe generation and acceleration sites. The distribution of 
Qnt for such high thermal fraction regions has a Gaussian 
shape, with mean Qnt ~ —0.78 and narrow spread of 0.16. 
It is interesting that this narrow distribution arises from 
multiple spatially resolved regions from an assorted set of 
galaxies with very different star formation histories. This is 
clearly indicative of a generic process of energy loss of CRe 
by synchrotron radiation and inverse Compton scattering 
as they propagate from their acceleration sites. The flatter 
part of Qnt (> —0.55) primarily arises from NGC 5236, in - 
dicating that ionization losses are dominant (|Longairll201ll ). 
iNiklas. Klein fc Wielebinskil (|l997l l. based on integrated fiux 
density of 74 galaxies found the average Qnt ~ —0.83 which 
is comparable to the Qnt value that we have obtained. The 
galaxy integrated spectral index is an average over tens of 
kpc where the flux density is highly dominated by bright 
regions which are mostly spiral arms or Hll regions. Our 
analysis shows that the contribution to the average spectral 
index primarily arises at least on kpc scales. 

The unflUed histogram, corresponding to low thermal 
fraction {flt^^ < 10%) regions, has a much larger spread 
ranging from Qnt ~ —0.4 to — 1.8. The steepening in the 
range —1.3 < Qnt < —1 arises from the interarm regions, 
and indicates energy losses through synchrotron emission 
and inverse Compton scattering. Qnt < —1.3, is mostly from 
low flux density regions, primarily in the outer parts of the 
galaxies, and may be affected due to missing flux density in 
the higher radio frequency maps. Note that we have aver- 
aged in a region 40", which for all these galaxies corresponds 
to ~ 1 kpc to 3 kpc, and hence the flatter part of Qnt distri- 
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bution may be affected due to bright regions. As mentioned 
earlier, the absolute value of Ont can change by ~10% due 
to uncertainty in flux scal e at 333 MHz. A recent study by 
iPaladino. Murgia fc Orrftl (flooa ') for NGC 628, NGC 3627 
and NGC 7331, using infrared emission (A70/^m) as a tracer 
of star forming regions, found that a obtained using VLA 
observations at 325 MHz and near 1.4 GHz, is flatter for 
infrared bright regions. This is similar to what we observe 
in our set of galaxies. The high angular resolution 333 MHz 
observations presented in this paper can be used to study 
in detail the spatially resolved FIR - radio and CO - FIR 
correlations, which will be reported in a future paper. 



APPENDIX A: THERMAL NONTHERMAL 
FLUX SEPARATION 



The lTabatabaei et all (|2007l ') method uses an Hq map of the 
galaxy as a template of thermal emission. Since Ha emission 
is affected by extinction, far-infrared data at A70/im and 
A160/im are used to obtain a model for extinction, which is 
then used to correct the Ha data. The corrected Ha image 
is then used to predict the thermal emission at radio fre- 
quencies of interest (the whole proces s is described in detail 
in sections 3 to 6 of lTabatabaei et al.li200 7V In this section, 
we summarize the various steps and the intermediate results 
that we obtained in the separation process. 

• Data sets: The far infrared images at A70^m and 
A160/im were obtained from the publicly available data 
from the Spitzer Inf rared Nearby Galaxy Survey (siNGS; 
iKennicutt et al.ll2003h . All A70pim images have a pixel size 
of 4.5" and a point spread function (PSF) of about 16". 
At A160/im, each pixel is 9" in size and has a PSF of 40". 
Both the X70^m and A160/im maps are calibrated in surface 
brightness units of MJy sr~^. 

The continuum subtracted Ha images for NGC 1097 
(1.5m CTIO, filter: CT6586) and NGC 6946 (KPNO 2m, 
filter: KP1563) were obtained from the ancillary data at the 
SINGS website. The maps were in units of DN s~^pixel~^, 
which was converted to erg s~^cm~'^ using the calibration 
provided in the SINGS Fifth Data Delivery documentatior|f|. 

The continuum subtracted Ha images for NGC 4736, 
observed with the I m JKT at La Palma, filter: Ha6570, 
l|Knapen et all 120041 '). NGC 5055 observed with the 2.3m 
telescope at KPNO, filter: 6580 and NGC 5236 observed 
with the 0.9m telescope at CTIO, filter:6563 were down- 
loaded from the NED. 

The counts per second (cps), in these maps were con- 
verted to apparent magnitude, mAB, using the zero- 
point given in the FITS header. The apparent magni- 
tude mAB was further converted to specific intensity, using 
/^(erg s-^cm-^Hz-i) = io-("ab+48.6)/2.6 rpj^g flux in units 
of erg s~^cm~^ was computed using, / — fvdu = fiyCj^. 

All the X70fim and continuum subtracted Hq maps were 
convolved to the coarser 40" PSF of A160/.im, and pixels 
were re-gridded to 9". All the subsequent calculations were 
done on a pixel by pixel basis. 



• Dust Temperature and Opt ical Depth: 

Following iTabatabaei et al.l (|2007l ). as a first step for ex- 
tinction correction of Hq maps we need to estimate the 
optical depth due to obscuring dust. For this the colour 
temperature of the dust (Tdust) is found by fitting a black 
body spectru m incorporating dust abso rption efficiency (see 
equation 1 of ITabatabaei et ahl |2007| ) to the far infrared 
maps at \70fim and A160/im. The distribution of the dust 
temperature for all the five galaxies is shown in the left 
panel of Figure lAll The mean temperature estimated for 
all the galaxies is between 19—22 K. Note that a sin- 
gle dust temperature model has been assumed. The op- 
tical depth at 16 0/im (rieoum) was then obtained by us- 
ing equation 2 of ITabatabaei et ahl (|2007l ) using the esti- 
mated Tdust. The Hq optical depth was obtained as rua ~ 
/d X 2200 X rim ,,^ llKriigell |2003| ) with Ha filling factor 
/d = 0.33 (Dic kinson et al.1 12OO3I '). For all the five galax- 
ies the THa lies in the range 0.03 - 0.6 at linear scales > 1 
kpc, with the highest values being towards the center of the 
galaxies. 

• Emission Measure: 

Extinction corrected Hq maps was used to estimate the 
emission measure [EM = Jn^ dl pc cm , where He is 
the thermal electron d ensity) usi ng equation 3 and 4 of 
ITabatabaei et al.l (120071 ) (also see IValls-Gabaudlll998l ) and 
assuming an electron temperature of Te — 10* K. The 
emission measure radial profile for all the five galaxies is 
shown in the middle panel of Figure IXT] The emission mea- 
sure towards the center of the galaxy lies in the range of 
- 10^ - 10" pc cm"*^ and falls off to ~ 10 - 10^ pc cm"*^ 
towards the edge of the galaxy. These range of emission mea- 
sure values a re consistent with those observed in the Milky 
way (see e.g. iBerkhuiisen. Mitra fc MuelleiibOOeD . 

• Thermal emiss ion at Radio frequen cies: 

Equation 5 and 6 of lTabatabaei et~aL iiiooj), were used to 

determine the radio continuum optical depth and the bright- 
ness temperature (Tb) using the Rayleigh- Jeans approxima- 
tion. The thermal flux density {Si,, th) at a radio frequency 
1/ was obtained from the Tb using, 



S,, A 



Jy beam^ 



8 



18xlO-(AfL) fll^) {-) 
Varcsec/ Varcsec/ VGHz/ 



Tb 
K 



Here, Sniaj and ^niin size of the restoring beam in the radio 
map (in this case 40"). 

For all the five galaxies, the pixel-wise distribution of the 
thermal fraction {fth = <5'i^,th/<S'i/,tot) at 333 MHz and near 
1 GHz are shown in the rightmost panel of Figure lAll The 
mean fth at 333 MHz and near 1 GHz for all the galaxies 
were found to be less than 5% and less than 12% respec- 
tively. However, in certain giant Hll regions, the thermal 
fracti on can go up to 10% at 333 MHz and 30% near 1 
GHz. ITabatabaei et ahl l|2007l ) pointed out that the primary 
source of error in the thermal fraction arises from the un- 
known value of the To. Using similar arguments, we estimate 
the uncertainty in the thermal fraction at 333 MHz and near 
1 GHz to be ~ 10% and 15% respectively. 
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Figure Al. Left panel: The pixel-wise distribution of the estimated dust temperature (T^ust)- The bins are at intervals of 0.5 K. Middle 
panel: The radial profile of the estimated emission measure {EM = / d/ pc cm ®). The EM is azimuthally averaged over annuli of 
one beamwidth, i.e, 40". Right panel: The pixel-wise distribution of the thermal fraction (/th) with a bin size of 0.5%. The solid black 
histograms show the distribution at 333 MHz, while the dashed red histograms show the distribution near 1.4 GHz. The top, middle and 
bottom rows are for NGC 1097, NGC 4736 and NGC 5055 respectively. Each pixel is of size 10", corresponding to a physical scale of 
about 700 pc, 400 pc and 340 pc respectively. 
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Figure Al. contd... Left panel: The pixel-wise distribution of the estimated dust temperature (Tdust)- The bins are at intervals of 0.5 
K. Middle panel: The radial profile of the estimated emission measure {EM = ^ nidi pc cm ^ ) . The EM is azimuthally averaged over 
annuli of one beamwidth, i.e, 40". Right panel: The pixel-wise distribution of the thermal fraction (/th) with a bin size of 0.5%. The 
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pc and 330 pc respectively. 



Baars J. W. M., Genzel R., Pauliny-Toth I. I. K., et al, 

1977, A&A, 61, 99 
Beck, R., Grave, R., 1982, A&A, 105, 192 
Beck, R., 1991, A&A, 251, 15 

Beck, R., Shoutenkov, V., Ehle, E., Harnett, J. I., Haynes, 
R. F., Shukurov, A., Sokoloff, D. D., Thierbach, M., 2002, 
A&A, 391, 83 

Beck, R., Fletcher, A., Shukurov, A., et al., 2005, A&A, 

444, 739 
Beck, R. 2007, A&A, 470, 539 
Bell, R., 1978, MNRAS, 182, 443 

Berkhuijsen, E. M., Beck, R., Hoernes, P., 2003, A&A, 398, 
937 

Berkhuijsen, E. M., Mitra, D., & Mueller, P. 2006, As- 

tronomische Nachrichten, 327, 82 
Biermann, P. L., Strom, R. G., 1993, A&A, 275, 659 
Bogdan, T. J., Volk, H. J., 1983, A&A, 122, 129 
Braun, R., Oosterloo, T. A., Morganti, R., Klein, U., Beck, 

R., 2007, A&A, 461, 455 
Briggs, D. S. 1995, Ph.D. Thesis, New Mexico Institute of 

Mining Technology, Socorro, New Mexico, USA 
Broadbent, A., Osborne, J. L., Haslam, C. G. T. 1989, 

MNRAS, 237, 381 
Brown, R. H., Hazard, C., 1961, MNRAS, 122, 479 
de Bruyn, A. G., 1977, A&A, 54, 491 



Carilh, C. L., Holdaway, M. A., Ho, P. T. P., de Free, C. 

G., 1992, ApJ, L59, 399 
Caswell, J. L., Wills, D., 1967, MNRAS, 135, 231 
Chyzy, K. T., Buta, R. J., 2008, ApJ, 677, L17 
Condon, J. J., 1987, ApJ SuppL, 65, 485 
Condon, J. J., 1992, ARA&A, 30, 575 
Condon, et al. 1998, AJ, 115, 1693 
Cornwell, T. J., Perley, R. A., 1992, A&A, 261, 353 
Dickinson, C., Davies, R. D., Davis, R. J., 2003, MNRAS, 

341, 369 

Duric, N., Dittmar, M. R., 1988, ApJ, 332, L67 
Gioia, L M., Gregorini, L., 1980, A&A SuppL, 41, 329 
Gioia, L M., Gregorini, L., Klein, U., 1982, A&A, 116, 164 
Green D.A., 1998, A Catalogue of Galactic Super- 
nova Remnants (1998 September version). MuUard 
Radio Astronomy Observatory, Cambridge, United 
Kingdom (available on the World-Wide- Web at 
http://www.mrao.cam.ac.uk/surveys/snrs/ ) 
Heesen, V., Beck, R., Krause, M., Dettmar, R. J., 2009, 

A&A, 494, 563 
Hummel, E., 1980, A&AS, 41, 15 
Hummel, E., Bosma, A., 1982, AJ, 87, 242 
de Jong, M. L., 1966, ApJ, 144, 153 
de Jong, M. L., 1967, ApJ, 150, 1 

Karachentsev, I. D., Sharina, M. E., Huchtmeier, W. K., 



18 Basu et al. 



2000, A&A, 362, 544 
Karachentsev, I. D., Sharina, M. E., Dolphin, A. E., et al., 

2002, A&A, 385, 21 

Karachentsev, I. D., Sharina, M. E., Dolphin, A. E., et al., 

2003, A&A, 398, 467 

Kellermann, K. I., PauUny-Toth, I. I. K., Williams, P. J. 

S., 1969, ApJ, 157, 1 
Kennicutt, Jr., R. C, Armus, L., Bendo, G., et al. 2003, 

PASP, 115, 928 
Klein, U., Emerson, D. T., 1981, A&A, 94, 29 
Klein, U., Beck, R., Buczilowski, U. R., Wielebinski, R., 

1982, A&A, 108, 176 
Klein, U., Grave, R., Wielebinski, R., 1983, A&A, 117, 332 
Klein, U., Wielebinski, R., Morsi, H. W., 1988, A&A, 190, 

41 

Knapen, J. H., Stedman, S., Bramich, D.-M., Folkes, S. L., 

Bradley, T. R., 2004, A&A, 426, 1135 
Knapik, J., Soida, M.,Dettmar, R.-J., Beck, R., Ur- 

banik, M., 2000, A&A, 362, 910 
Kothes, R., Fedotov, K., Foster, T. J., Uyaniker, B., 2006, 

A&A, 457, 1081 
Kriigel, E. 2003, The physics of interstellar dust. The 

physics of interstellar dust, ed. E. Kruegel (loP Series in 

astronomy and astrophysics, Bristol, UK: The Institute of 

Physics). 

Kuril'chik, V. N., Onishchenko, L. V., Turyvskii, V. M., 
1967, Sov. Astron. J., 11, 528 

Kuril'chik, V. N., Andrievskii, A. E., Ivanov, V. N., Span- 

genberg, E. E., 1970, Sov. Astron. J., 13, 881 
van der Kruit, P. C., Allen, R. J., Rots, A. H., 1977, A&A, 

55, 421 

Longair, M. S., 2011, High Energy Astrophysics 3-ed, Cam- 
bridge University Press. 

Maddox, L. A., Cowan, J. J., Kilgard, R. E., Lacey, C. K., 
Prestwich, A. H., Stockdale, C. J., Wolfing, E., 2006, AJ, 

132, 310 

Mayya, Y. D., Carrasco, L., Luna, A., 2005, ApJ, 628, L33 

Murgia, M., Heifer, T. T., Ekers, R., Blitz, L., MoscadeUi, 
L., Wong, T., Paladino, R., 2005, A&A, 437, 389 

Murphy, E. J., Helou, G., Condon, J. J., Schinnerer, E., 
Turner, J. L., Beck, R., Mason, B. S., Chary, R.-R., Ar- 
mus, L., 2010, ApJ, 709, L108 

Nagar, N. M., Falcke, H., Wilson, A. S., 2005, A&A, 435, 
521 

Niklas, S., Beck, R., 1997, A&A, 320, 54 
Niklas, S., Klein, U., Wielebinski, R., 1997, A&A, 322, 19 
Ondrechen, M. P., van der Hulst, J. M., 1983, ApJ, 269, 
L47 

Ondrechen, M. P., van der Hulst, J. M., Hummel, E., 1985, 

AJ, 90, 1474 
Ondrechen, M. P., 1989, ApJ, 342, 39 
Paladino, R., Murgia, M., Orru, E., 2009, A&A, 503, 747 
Rieke, G., et al., 2004, ApJS, 154, 25 
Roy, S., Rao P. A., 2004, MNRAS, 349, L25 
Schommcr, R. A., Sullivan, W. T., 1976, Ap. Letters, 17, 

191 

Seaquist, E. R., Odegard, N., 1991, ApJ, 369, 320 
Segalovitz, A., 1977a, A&A, 61, 59 

Segalovitz, A., 1977b, A&A, 54, 703 

Smith, B. J., Lester, D. F., Harvey, P. M., Pogge, R. W., 

1991, ApJ, 380, 677 
Steer, D. G., Dewdney, P. E., Ito, M. R., 1984, A&A, 137, 



159 

Sukumar, S., Klein, U., Grave, R., 1987, A&A, 184, 71 
Swarup, G., Ananthakrishnan, S., Kapahi, V. K., et al. 

1991, Current Science, 60, 95 
Tabatabaei, F., Beck, R., Krugel, E., , Krause, M., Berkhui- 

jsen, E. M., Gordon, K. D., Menten, K. M. 2007, A&A, 

475, 133 

Valls-Gabaud, D., 1998, PAS A, 15, 111 

de Vaucouleurs, G., et al. 1991, Third Reference Catalogue 

of Bright Galaxies (Berlin: Cambridge University Press) 
Wong, T., Blitz, L., 2001, ApJ, 540, 771 
Wright, A. E., Wark, R. M., Troup, E., Otrupcek, R., Hunt, 

A., Cooke, D. J., 1990, PASA, 8, 261. 



